arXiv:1505.01814vl [astro-ph.SR] 7 May 2015 


Mon. Not. R. Astron. Soc. 000, ^|9] (2015) Printed 8 May 2015 (MN style file v2.2) 


The Central Star Candidate of the Planetary Nebula 
Sh2-71: Photometric and Spectroscopic Variability 

T. Mocnik, 1,2 * M. Lloyd , 3 D. Pollacco 4 and R. A. Street 5 

1 Isaac Newton Group of Telescopes, Apartado de Correos 368, E-38700 Santa Cruz de La Palma, Spain 

2 Astrophysics Group, Keele University, Staffordshire ST5 5BG, UK 

3 Jo dr el l Bank Centre for Astrophysics, School of Physics and Astronomy, University of Manchester, Manchester M13 9PL, UK 

4 Department of Physics, University of Warwick, Coventry CVf 7AL, UK 

5 Las Cumbres Observatory Global Telescope Network, 67f() Cortona Drive, Suite 102, Goleta, CA 93117, USA 


Accepted xxxx xxxxxxxx xx. Received xxxx xxxxxxxx xx; in original form xxxx xxxxxxxx xx 


ABSTRACT 

We present the analysis of several newly obtained and archived photometric and 
spectroscopic datasets of the intriguing and yet poorly understood 13.5-mag central 
star candidate of the bipolar planetary nebula Sh2-71. Photometric observations con¬ 
firmed the previously determined quasi-sinusoidal lightcurve with a period of 68 days 
and also indicated periodic sharp brightness dips, possibly eclipses, with a period of 
17.2 days. In addition, the comparison between U and V liglitcurves revealed that 
the 68-day brightness variations are accompanied by a variable reddening effect of 
A E(U — V) = 0.38. Spectroscopic datasets demonstrated pronounced variations in 
spectral profiles of Balmer, helium and singly ionised metal lines and indicated that 
these variations occur on a time-scale of a few days. The most accurate verification to 
date revealed that spectral variability is not correlated with the 68-day brightness vari¬ 
ations. The mean radial velocity of the observed star was measured to be ~26 km s' 1 
with an amplitude of ±40 km s' 1 . The spectral type was determined to be B8V through 
spectral comparison with synthetic and standard spectra. The newly proposed model 
for the central star candidate is a Be binary with a misaligned precessing disc. 

Key words: planetary nebula: individual: Sh2-71 - stars: variables: general - stars: 
emission-line, Be - binaries: close. 


1 INTRODUCTION 

Planetary nebula (PN) Sh2-71 (PN G035.9-01.1, a = 
19 h 02 m 00.3°, S = +2°09Tl" ,12000), see Fig. 1, was dis- 


nebula. The largest photometric study to date was presented 


covered by Minkowski (1946) who classified it as a diffuse 


and peculiar nebulosity. In 1959 the second Sharpless cata¬ 
logue of emission nebulae listed the object as a possible PN 
( Sharpless|1 959). Modern catalogues classify the object as a 
true PN ( Kohoutek|2001 1. 

The first extensive photometric observations of the vari¬ 
able 13.5-mag apparent central star of Sh2-71 were per¬ 
formed by |Kohoutek| between 1977 and 1979. |Kohoutek| 
(1979) reported brightness variability of more than 0.7 mag 


with indications of a smooth sine-like lightcurve. He also 
used the star’s colour indices to initially classify the star as a 
B8V and noted that the spectral type B8V cannot be respon¬ 
sible for the observed high excitation in the surrounding PN. 
Therefore, he suggested that the central star is a binary sys¬ 
tem in which the secondary star is hot enough to ionize the 
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by Mikulasek et al. (20071. They obtained 2,004 photomet¬ 
ric measurements in V, Rc and Ic band, and combined 
them with existing photometric data from |Kohoutek| ( |1979| ) 
and Jurcsik (19931. From the combined 4,268 measurements, 


spanning over 28 years, they identified the brightness vari¬ 
ations to be a quasi-sinusoidal lightcurve with a period of 
68.101 ± 0.010 days, and that the lightcurve minima flatten 


occasionally. Mikulasek et al. (2007) also noted pronounced 


spectral variations between two spectra, taken at different 
brightness phases. 

A detailed model for the central star as the progenitor 
of the collimated bipolar nebula was suggested by |Cuesta| 
& Phillips (1993). Primarily based on their four slit spectra 


covering both star and nebula and taken on two consecu¬ 
tive nights they suggested that the central star is a close 
binary system in which the material is being transferred 
from the main-sequence star to an ionising star. They fur¬ 
ther speculated that the material expelled through the sec¬ 
ond Lagrangian point creates the circumstellar disc, which 
collimates the fast stellar wind into a bipolar shape. Based 
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Figure 1 . RGB image of planetary nebula Sh2-71 and zoomed central region of the nebula (top right panel), identifying central star 
candidates CS1 and CS2. Note the nebular feature which is extending eastwards from CS1. The two comparison stars, which were used 
for the differential photometry, are labelled 1 and 2. Ho: (red), Sloan r’ (green), [O ill] (blue). 


on synthetic composite spectra, Cuesta & Phillips (19931 
determined the spectral type of the cooler, brighter compo¬ 
nent to F0V-A7V and estimated the effective temperature 
of the ionising star to be 10 5 K. To date, there is no direct 
evidence for a hot, ionising binary component. 

Although the variable 13.5-mag star has historically 


been identified as the central star (star CS1 in Fig. 1), Frew 
& Parker] ( |2007| ), have suggested however that a nearby, 
much dimmer and bluer, 19-mag star may be the true cen¬ 
tral star (star CS2 in Fig. l). |De Marco e t ah (2008) have 
concluded that there is insufficient evidence either way to 
give a definitive identification of the true central star. 

In this paper we present our imaging of the nebula and 
challenge the alternative suggestion for the central star. Fur¬ 
thermore, we present our analysis of new photometric and 
spectroscopic observations, coupled with some archived and 
yet unpublished spectroscopic datasets. From these obser¬ 
vations we derive several new parameters of the brighter 
central star candidate and refine certain previously deter¬ 
mined properties. We suggest a new model for the central 
star candidate in order to establish a link between the ob¬ 
served properties and the collimated bipolar shape of the 
PN. 


2 OBSERVATIONS AND DATA REDUCTION 

We have imaged Sh2-71 in broad- (Sloan r’, 3x30 s) and two 
different narrowband filters (Ha and [O ill], both 3x120 s) 
on 2013 October 21 with the Wide Field Camera (WFC), 


which is mounted on the 2.5-m Isaac Newton Telescope 
(INT) on La Palma in the Canary Islands. The images were 
processed in the usual way using theli, fits liberator and 
matlab to produce the RGB image shown in Fig. 1. 

The three photometric datasets of the ~ 13.5-mag cen¬ 
tral star candidate CS1, presented in this paper, were ob¬ 
tained with the 0.8-m telescope at Byrne Observatory at 
Sedgwick (BOS) in California, the 2-m Faulkes Telescope 
South (FTS) in Siding Spring, Australia, and FTS’s twin 
telescope, the 2-m Liverpool Telescope (LT) on La Palma. 
Observational details are listed in Table 1. All three pho¬ 
tometric datasets have been bias subtracted and flat field 
corrected by the automatic pipelines, provided by the in¬ 
dividual observatories. Combining the LT images and dif¬ 
ferential aperture photometry were carried out using the 
Starlink application packages ccdpack and photom. The 
non-variable comparison star, which was present in every 
image and yielded the highest signal to noise ratios of the 
measured differential magnitudes for CS1 in FTS and LT 
datasets, is shown in Fig. 1 as comparison star number 1. 
Because the BOS images were obtained with much longer 
exposure times compared to FTS and LT, the comparison 
star number 1 lay outside the linear range of the CCD’s full 
well capacity for some of the BOS images. Thus, the BOS 
differential photometry was measured with a fainter com¬ 
parison star which is labelled as comparison star number 2 
in Fig. 1. The differential photometry was performed using 
circular photometric apertures and annulus sky subtraction. 
The aperture radii were set to be equal to one average stellar 
full width at half maximum (FWHM) as measured for a set 
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of five most suitable stars in each image. Sky annulus was 
placed between two and three FWHM from the stellar cen¬ 
tres. This photometric configuration was found to maximize 
the signal to noise ratios of the differential magnitudes. We 
were only able to measure the differential photometry for 
CS1; CS2 was too faint to be detectable in the FTS and LT 
datasets and below the 5-sigma detection limit for the ma¬ 
jority of the images in the BOS dataset. The log of the final 
differential magnitudes for CS1 is given in Online Appendix 
in Tables 1.1-1.5. 

The first spectroscopic dataset comprises 25 unpub¬ 
lished, archived slit spectra of the central star candidate CS1 
which were obtained from the Intermediate Dispersion Spec¬ 
trograph (IDS), mounted on the INT. A further 9 archived 
and unpublished high-cadence and high-resolution echcllc 
spectra were obtained from the University College London 
Echelle Spectrograph (UCLES) on 3.9-m Anglo-Australian 
Telescope (AAT) in Siding Spring, Australia. The third 
spectroscopic dataset, purposely obtained for this research 
project, was obtained from the Fibre-Fed Robotic Dual- 
Beam Optical Spectrograph (FRODOSpec) at the LT. Ob¬ 
servational details for the individual spectroscopic datasets 
are presented in Table 2. Pre-processing, extraction of spec¬ 
tra and wavelength calibration of the INT dataset were 
performed with STARLINK application packages FIGARO and 
echomop. The AAT spectra were reduced with the iraf 
application package echelle. The LT spectra were reduced 
with the fully automated data reduction pipeline provided 


by Barnsley et al. (20121. 


3 RESULTS 
3.1 Imaging 

Our RGB image of Sh2-71 is shown in Fig. 1. As noted in 
the Introduction, the bright star CS1, historically believed 
to be the central star of Sh2-71, is not hot enough to ionise 
the nebula. This led Frew fe Parker| |2007;) to propose that 
the faint (19 th mag) blue star, offset 7 arcsec NW of CS1 and 
identified here as CS2, may be the true central star. They 
claim that the measured brightness of CS2 corresponds to 
the expected luminosity of the progenitor’s remnant at the 


nebula’s distance of 1 kpc. In addition, Frew & Parker (2007) 


note that CS2 is also closer to the apparent geometric centre 
of the nebula, whereas CS1 is slightly offset from the centre. 
The WFC images of the nebula have been used to investigate 
the nebular features close to the centre of the nebula to try 
to find evidence as to which star is the true central star of 
the PN. 

The central position of CS2 favours its identification 
as the central star, although several other PNe have been 


discovered to possess off-centre central stars (De Marco 


2009). However, the somewhat symmetrical nebular struc¬ 
ture about a line eastwards from CS1 may suggest a physical 
link to this star. Furthermore, since the most likely colli¬ 
mating mechanism for the bipolar shape of a PN is a binary 


central star (De Marco 2009), the progenitor star of the bipo¬ 


lar PN Sh2-71 is expected to have a companion star which 
would provide additional brightness contribution. Thus, a 
brighter central star candidate, CS1, seems a more plausi¬ 
ble choice for a binary central star hypothesis. In such a 


case, CS1 would comprise a hot ionising component and a 
cooler, optically brighter companion which would dominate 
the observed optical flux, effectively hiding the hot ionising 
companion. A similar binary system has been identified as 
the central stellar system of another bipolar PN, NGC 2346 


(Mendez & Niemela (19811, also see subsection 4.2). Finally, 
the reported mean RV of the PN and measured mean stellar 
RV of CS1 are very similar (see subsection 3.3.3), as would 
be expected for a central star and its planetary nebula, al¬ 
though as yet there is no RV measurements for CS2. Thus 
CS1 remains a reasonable candidate for the central star at 
this stage. 

3.2 Photometry 

The complete detailed photometric analysis is presented in 


Mocnik (2013). In this paper the main photometric results 


are summarised in Fig. 2, which shows a set of four individ¬ 
ual lightcurves obtained from three different telescopes. The 
BOS B, FTS V and LT V datasets are used to reveal the 
shape of the lightcurve, whereas the overlay between LT V 
and LT U illustrates the difference in axes scaling between 
differential magnitudes AV and A U, which can be used to 
estimate the reddening effect. Inconsistencies in the filters 
used and the available comparison stars for each image make 
it difficult to compare these photometric measurements di¬ 
rectly. Here, the BOS data have been obtained using com¬ 
parison star 2 (see Fig. 1) and the FTS and LT data have 
been obtained using comparison star 1. Thus the FTS and 
LT datasets are directly comparable, however, the BOS dif¬ 
ferential magnitudes had to be offset. Fig. 2 also shows the 


predicted lightcurve given by Mikulasek et al. (2007). 

First of all, we can see from Fig. 2 that the com¬ 
bined lightcurve is in rough agreement with the period 


and phase reported by Mikulasek et al. (20071 (JD max = 


2 449 862.02±0.32+68.101±0.010(£-96)). The comparison 
confirms that the photometric nature of CS1 has remained 
the same since the last Mikulasek’s observation in 2002. The 
disagreement between the observed and predicted decreas¬ 
ing part of the lightcurve is most likely due to the current 
deviation of the lightcuve’s shape from the assumed sinu¬ 
soidal brightness variations. 

Secondly, several sharp dips can be seen in the 
lightcurve. These are significantly larger than the uncertain¬ 
ties of individual data points but much smaller than the al¬ 
most 1-mag peak-to-peak variation in the overall lightcurve. 
Attempts at a formal periodic analysis (Mocnik 2013) were 
unsuccessful due to insufficient data. However, a simple 
visual fit of these dips, made by trying various periods 
and epochs, revealed that the five most pronounced bright¬ 
ness dips can be characterised as periodic with a period of 
17.2±0.1 days. The overlap between the brightness dips and 
best-fit epoch and period is illustrated with the periodic ver¬ 
tical lines in Fig. 2. Such a potentially periodic behaviour 
points to the intriguing possibility of an eclipsing binary 
system, however, many more data points would be required 
for a more reliable and accurate analysis of these possibly 
periodic brightness dips. 

Lastly, the difference between the axes scales for LT V 
and U filter in Fig. 2 reveals that the brightness variation 
amplitude is larger for U filter and therefore that the cen¬ 
tral star exhibits a variable reddening effect along with the 
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Table 1. Observational details for photometry. The table lists the names of the 
telescopes, the dates of the first and last observation, applied filters, exposure times, 
mean signal to noise ratios for the ~ 13.5-magnitude CS1 and the final number of 
useful photometric data points. 


Telescope 

Duration of 
observations 

Filter 

Exposure 
time (s) 

S/N 

Number of 
data points 

BOS 

2011 Jun 26 - 2011 Jul 28 

B 

600 

543 

21 

FTS 

2011 Oct 4 - 2012 Mar 23 

B 

15 

92 

7 



V 

15 

157 

7 

LT 

2012 Mar 31 - 2012 Jul 20 

U 

3x24 

64 

34 



V 

3x3 

171 

36 


Table 2. Overview of spectroscopic observations. The INT exposure time is given as a median value. Mean 
signal to noise ratios are calculated for central wavelengths, except for AAT, where it is provided for H/3 region. 


Telescope 

Duration of 
observations 

Wavelength 
range (A) 

Resolving 

power 

Exposure 
time (s) 

S/N 

Number of 
spectra 

INT (IDS) 

1996 Jul 10 - 1999 May 4 

3820-4800 

4500 

1500 

161 

25 

AAT (UCLES) 

2000 May 21 

3700-5100 

44000 

1800 

20 

9 

LT (FRODOSpec) 

2012 Jun 16 - 2012 Jul 19 

3900-5100 

5500 

1200 

16 

2 



5900-8000 

5300 

1200 

31 

3 



Figure 2. The comparison between predicted and measured lightcurves of CS1 with tentative indication of sharp periodic brightness 
dips. For compatibility reasons, the BOS differential magnitudes were offset by -3.4 mag. The dashed line represents the predicted 
lightcurve by |Mikulasek et al.| < |2007| . Potentially periodic brightness dips are indicated by engraved vertical lines with a period of 17.2 
days. The difference in axes scaling between LT V (black) and U (blue) lightcurves suggests a variable reddening effect. The data points 
are connected with lines for easier presentation purposes. 
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Wavelength (Angstrom) 


Figure 3. LT spectra of CS1 near the H/3 profile. The two spec¬ 
tra demonstrate the pronounced variations of H/3 (4861 A) and 
Fell lines (4924 and 5018 A). The first spectrum, plotted with 
the black colour, has been obtained on 2012 July 1 with a corre¬ 
sponding 68-day brightness phase of 0.754. The second spectrum, 
shown in red, was obtained on 2012 July 19 with a phase of 0.018. 

brightness variations. The variable component of reddening 
effect was estimated to be AE(U — V) = 0.38. On the other 
hand, AE(B — V) could not be estimated reliably because 
BOS B lightcurve did not chronologically overlap with any 
other photometric dataset, and because of insufficient num¬ 
ber of data points in FTS B dataset. 


3.3 Spectroscopy 

3.3.1 Spectroscopic Variability 



Note the absence of pronounced spectroscopic variations in the 
4.5 hour time-span and that H/3 emission components are reversed 
compared to the LT spectra in Fig. 3. 

spectral variations in the LT and INT datasets and their ab¬ 
sence in the AAT dataset imply that the spectral variations 
must occur on a time-scale of a few days. 

3.3.2 Spectral Type 


Fig. 3 shows two fully reduced and normalized spectra of 
the observed star near the H/3 line, obtained with the LT at 
different phases of the 68-day lightcurve. Pronounced spec¬ 
tral variations of the emission components of Balmer and 
singly ionised metal line profiles confirm the variations, re¬ 
ported byjMikulasek et al. ( |2007[ ). However, in contrast to 
Mikulasek’s reporting we did not notice pronounced [O ill] 
lines in any of the retrieved spectral datasets. Possibly the 
spectra analysed by |Mikulasek et al.| were contaminated by 
the nebular [O ill] line. 

We employed a dataset of 25 archived INT slit spec¬ 
tra to provide the most accurate verification to date if the 
spectral variations are correlated to 68-days brightness vari¬ 
ations. For each of the reduced and normalized spectrum 
we have calculated the corresponding brightness phase from 


We compared the measured spectra of the central star can¬ 
didate CS1 with synthetic spectra, computed by Munari; 


the ephemeris, published by Mikulasek et al. (2007). Nor¬ 


malized INT spectra, arranged according to their calculated 
brightness phase, are shown in Fig. 4. 

In contrary to expressed expectations in the previous 
papers related to Sh2-71, Fig. 4 demonstrates that the 
spectral variations are not correlated with 68-days bright¬ 
ness variations. Due to a large chronological spread of the 
archived INT spectra we were unable to determine any other 
potential period in the spectral variability. For this task the 
spectra would have to be obtained with a higher cadence, 
preferably with a daily cadence. On the other hand, 9 high- 
resolution echelle spectra from the AAT with a cadence of 30 
minutes did not reveal any noticeable spectral variations in 
the time-span of 4.5 hours (see Fig. 5). The presence of the 


et al. (2005). The best-fitted stellar atmospheric param¬ 
eters were determined to be T e g 11000-13000 K, log(<?) 
4.0 ± 0.25 cm s’ 2 , i> ro t • sin(i) 200^25 km s’ 1 with an indi¬ 
cated high value of metallicity. These parameters classify the 
observed star as a B9-B7 main-sequence star. The further 
comparison to standard spectra, provided by |Bagnulo et ah] 
12003), confirmed a late B main-sequence spectral type. The 
comparisons and details regarding the synthetic and stan¬ 
dard spectra are presented in Online Appendix and |Mocnik| 
( 20131 ). 

Our ~B8V spectral type classification of CS1 is in 
agreement with earlier B8V classification by |Kohoutek| 
(1979), who determined the spectral type by colour indices. 
On the other hand, our ~B8V classification refutes the spec¬ 


tral type A7V-F0V by Cuesta & Phillips (19931 who were us¬ 
ing slit spectra in the wavelength range 6000-7400 A. They 
noted, however, the presence of Si II lines (6347 and 6371 A) 
in their spectra, which correspond to stars with tempera¬ 
tures higher than A7V. The spectral wavelength range, cov¬ 
ered by this project’s INT spectra, has revealed many more 
spectral features that cannot be explained with an A7V 
star, such as the presence of Sill lines (3856 and 3863 A), 
He I (4026 and 4471 A) and absence of Mm (4033 A), Cal 
(4227 A) and Till (4468 A). The synthetic spectra fitting 
and comparison with the standard spectra provide strong 
evidence that the spectral type of the central star candidate 
CS1 is ~B8V. 
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Figure 4. A set of INT spectra, ordered by their phase within the 68-day photometric variability. The phases of individual observing 
nights are provided under each corresponding date. The most pronounced spectral variability may be noticed for Balmer series, singly 
ionised metal lines and neutral helium (4026 and 4471 A). Note that the spectral variability is not correlated with the brightness phase. 


© 2015 RAS, MNRAS 000, [TJ^ 









The Central Star Candidate of the Planetary Nebula Sh2-71 7 


3.3.3 Radial Velocity 

Radial velocities (RVs) of CS1 in the INT and AAT datasets 
have been measured in two steps. First, we measured the 
relative RVs with respect to the first spectrum in each of 
the two datasets. For the INT dataset we used the star- 
link tool HCROSS, the template INT spectrum of CS1 from 
1996 Jul 10 and a cross-correlation wavelength range 3960- 
4600 A. For the AAT dataset we used the iraf tool fx- 
COR, template AAT spectrum of CS1 with mid-exposure 
time 15:22:20 UT and a cross-correlation wavelength range 
4000-5000 A. Additional binning was not applied to either 
of the two datasets. The second step was to convert relative 
RVs into absolute values. This was done for both datasets by 
measuring the wavelength of Mgn 4481 A, which has been 
identified as the most stable and suitable individual spectral 
line for measuring RVs in our spectra. The final absolute he¬ 
liocentric RVs for the INT and AAT datasets are listed in 
Online Appendix in Tables 2.1 and 2.2, respectively. 

Our INT RV measurements revealed that the mean RV 
of the central star candidate CS1 is ~26 km s" 1 and that 
RV varies with an amplitude of ±40 km s" 1 (see top panel 
in Fig. 6). Due to the large spread in observation time of 
our INT spectra, we were unable to determine the period 
of a potential periodicity of the RV variations (see bottom 
panel in Fig. 6). However, the cadence and the number of 
data points were sufficient to reject any correlation with the 
68-days brightness variations (see middle panel in Fig. 6). 
In addition, because of the absence of RV variations in AAT 
spectra with 30-minutes cadence (see Table 2.2 in Online 
Appendix) we also rejected hourly time-scales for a potential 
RV periodicity. 

As discussed in subsection 3.1, some researchers have 
suggested that CS1, the star observed and presented in this 
paper, is actually not the central star of the PN (see Fig. 1). 
However, a useful result, derived from our RV analysis, is 


agreement with the reported mean RV of the PN, 24.7 km s" 


I Schneider et al. 1983), as would be expected for a central 


star and its surrounding PN. We note that the RV of the 
much dimmer 19-magnitude CS2, an alternative candidate 
for the central star, is as yet unknown. 


4 SUGGESTED MODEL FOR THE CENTRAL 
STAR 

4.1 Cataclysmic binary 


The cataclysmic binary model, suggested by [Cuesta ~Ez\ 
Phillips (19931, provides a plausible explanation for the ob¬ 


served Ha spectral line profiles as well as for the bipolar 
shape of the PN. However, a small separation between the 
stars, required for mass transfer to establish between a B8V 
and an ionising star, would correspond to an orbital period 
of less than 1.0 day. This would cause the RV of Balmer 
emission components to vary with ± >250-sin(i) km s. 
Our set of AAT spectra did not reveal any noticeable spec¬ 
tral line profile or RV variations of Balmer emission lines 
in the 4.5-hour time-span which demonstrates that the cat¬ 
aclysmic binary model is highly unlikely. In addition, this 
model does not provide a plausible explanation for the 68- 
days brightness variation period. 



Figure 6. INT RV measurements. Top panel: RV versus time. 
Dashed grey line indicates the mean RV of 26.3 km s" 1 . Middle 
panel: RV versus 68-day brightness phase. Note the absence of any 
correlation between the RV variations and the 68-day brightness 
phase. Bottom panel: Lomb-Scargle periodogram. Note that the 
periodogram does not reveal any RV periodicity in the frequency 
range from 0 and 1 day -1 . 


Be binary with precessing disc 


Our analysis of several potential models revealed that the 
observed spectrophotometric properties of CS1 are best re¬ 
produced with a close Be binary system with a misaligned 
precessing circumstellar disc. Other models such as Herbig 
Be stars, magnetically warped circumstellar discs, stellar 


pulsations etc. are discussed in Mocnik (2013). 


First of all, near edge-on rotating circumstellar discs 
around Be stars produce profiles of Balmer and singly 
ionised metal spectral lines similar to ones observed in our 


spectra. Hanuschik (19891 provided the relation between 


the full width at half maximum (FWHM) of the Balmer 
lines and the stellar rotational velocity. Our observed mean 
FWHM of Ha (9.3 A) and H/3 (5.7 A) emission profiles yield 
a stellar rotational velocity of ~270 km s" 1 , which is in rough 
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agreement with the measured rotational velocity of the cen¬ 
tral star (200^25 km s_1 > see subsection 3.3.2). 

Secondly, the observed spectral variability can be ex¬ 
plained by various contributing factors, which have been re¬ 
ported for a significant fraction of Be stars, e.g. variable wind 
absorption, nonradial stellar pulsations, starspots, magneti¬ 


cally modulated circumstellar rotation, etc. (Porter & Riv- 


(20031, Rivinius et al. (20131). In a close Be binary 
system we might also expect to see the periodic spectral and 
RV variations arising from the orbital motion of the ionising 
star and consequent periodic heating of outer Be disc layers, 
similar to tp Per, which is a binary between B0.5IVe and an 
O-type subdwarf (sd06) ( |Bozic et al.|1995 l. Poeckert (19811 
reported that spectra of <p Per also exhibit periodic inter¬ 
changing blue- and red-dominated Balmer emission compo¬ 
nents, as indicated by Figs. 3 and 5 in the case of CS1. We 
do not see evidence that the spectral variations are peri¬ 
odic, however, our data were only sensitive to very short 
timescales (<1 day) and were sufficient to reject also the 
68-day periodicity. 

The expected orbital period of a close Be binary would 
be longer than in the case of a cataclysmic binary model. 
Having only two reasonable assumptions available, i.e. that 
the primary star had a mass higher than the B8V companion 
star in the pre-PN phase, ~2.6 M@ (Belikov[l995 I, and that 
the Roche lobe was filled during the red giant phase of the 
primary star, we can only derive rough theoretical lower and 
upper limits to the current orbital period to be 1 and 300 
days, respectively. However, a similar close binary system of 
an early main sequence star (A5V) and an O-type subdwarf 
was found in the centre of the bipolar PN NGC 2346 with an 
orbital period 16 days ( Mendez &; Niemela|198l] ). If a similar 
orbital period is applied to Sh2-71, this would comply with 
the absence of spectroscopic variations in hourly time-scales 
of AAT dataset and presence of variations in daily time- 
scales in the INT and LT datasets. A similar orbital period 
would also coincide with the tentatively identified eclipsing 
period of 17.2 days. 

Furthermore, the mass transfer in the pre-PN evolu¬ 
tion, required for the formation of a subdwarf, provides an 
explanation for the observed high metallicity in the spec¬ 
tra of the observed star, because the metal-rich material, 
recently dumped onto the B8V star in the pre-PN phase 
would not have had enough time to be fully mixed into the 
main-sequence star (Smal ley|1997 ). 

The photometric behaviour of the observed star can be 
explained by a variety of different contributing factors, such 
as starspots, stellar pulsations, reflection and/or irradiation 
if in a close binary system etc., but for the vast majority of 
cases, the effect would be insufficient or it would be corre¬ 
lated with spectral variations. On the other hand, a close 
Be binary system with a precessing Be disc has been iden¬ 
tified as a possible plausible mechanism for the observed 
photometric properties and uncorrelated spectral variations 
of CS1. If the circumstellar Be disc is inclined with re¬ 
spect to orbital plane, the close binary ionising star can 
invoke tidal torques and cause the inclined disc to precess. 
The resulting periodic occultations of the companion star 
by the precessing disc can then be responsible for the ob¬ 
served 0.8-magnitude quasi-sinusoidal brightness variations. 
Furthermore, the detected variation of the reddening effect, 
A E(U — V) = 0.38, which coincides with the period and 


phase of the CSl’s lightcurve, corresponds to an estimated 
absorption of 0.7-0.8 mag in V-band and thus confirms that 
the extinction due to the precessing Be disc could be a real¬ 
istic explanation for the observed photometric behaviour. 
The quasi-sinusoidal lightcurve with a period of 68 days 
would then correspond to a precession of the circumstel¬ 
lar disc, whereas the uncorrelated spectroscopic variations 
would result from the various intrinsic contributing factors 
listed above or, if periodic, possibly a result from orbital 
motion of the binary components. 

The reported fraction of B8e among all B8 stars is 10 
per cent, the vast majority of them being main sequence 
stars (Kogure et al. 1982), whereas the binary fraction of Be 
stars is ~l/3 ( |Porter fc Rivini us 2003). Combined fractions 
suggest that about 3 per cent of B8 stars are B8Ve binaries. 
From these, Be+sdO binaries in particular are supposed to 
exist in abundance, yet because they are not easily distinct, 


only three such systems are known with certainty (Rivinius 
|et al.|[20 13). This provides a small but still realistic prob¬ 
ability that the central star candidate CS1 is a B8Ve+sdO 
binary. 

The proposed precessing Be-disc model is highly tenta¬ 
tive, however, it provides the explanation for the observed 
brightness variations and uncorrelated pronounced spectral 
variability. In addition, it also provides a potential collima- 
tion mechanism, required for the bipolar shape of Sh2-71. 
An extended daily-cadence sample of Ha and H/I profiles 
would be required to verify any potential periodicity of the 
spectral variations and thus provide further evidence for the 
Be binary hypothesis. 


5 CONCLUSIONS 

Using photometric datasets from the central star candidate 
CS1 of the planetary nebula Sh2-71 from the Byrne Ob¬ 
servatory at Sedgwick (BOS), the Faulkes Telescope South 
(FTS) and the Liverpool Telescope (LT), we tentatively 
identify the presence of periodic sharp brightness dips, pos¬ 
sibly eclipses, with a period of 17.2 ±0.1 days. We have also 
measured a periodic variable reddening effect with an am¬ 
plitude A E(U — V) = 0.38, which coincides with the 68-days 
period and phase of the star’s quasi-sinusoidal lightcurve. 

Spectroscopic datasets from the LT, the Isaac New¬ 
ton Telescope (INT) and the Anglo-Australian Telescope 
(AAT) confirmed pronounced variations in the shapes of the 
Balmer, helium and singly ionised metal lines. Our analy¬ 
sis revealed that these spectral variations are not correlated 
with the 68-days brightness variations and that the spectral 
variations occur on the time-scale of a few days. Compar¬ 
ing the observed with synthetic spectra yielded T e g 11000- 
13000 K, log (g) 4.0 ± 0.25 cm s' 2 , v TO t • sin(i) 200^25 km s" 1 
and indicated a high value of metallicity. The fitted parame¬ 
ters and the comparison with standard spectra classified the 
observed star as ~B8V. The radial velocity (RV) was mea¬ 
sured to vary with an amplitude of ±40 km s' 1 around the 
mean value of ~26 km s' 1 . The near agreement between the 
mean RV of the observed star and the reported mean RV of 
the nebula provides supporting evidence that the brighter 
candidate for the central star, CS1, indeed could be the true 
central star of Sh2-71. In addition, imaging with the INT 
revealed nebular features extending eastwards from CS1, 
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which indicates that the observed star is likely to be physi¬ 
cally linked with the planetary nebula. 

We have revisited and rejected the cataclysmic model 
for the central star, proposed by Cuesta & Phillips (19931. 
The newly obtained spectrophotometric properties of the 
observed star were best reproduced by a close Be binary 
system with a precessing circumstellar disc. However, mainly 
due to the lack of spectra with a daily cadence, this newly 
suggested model should be regarded as tentative. 
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The Central Star Candidate of the Planetary Nebula Sh2-71: 
Photometric and Spectroscopic Variability - Online Appendix 
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Abstract 

This online appendix lists the differential magnitudes and radial velocity measurements of the central 
star candidate CS1 and expands the spectroscopic analysis, introduced in the main article. Additional data, 
explanation and figures are focused in four sections: Differential Photometry, Radial Velocities, Synthetic 
Spectra and Standard Spectra. 


1 DIFFERENTIAL PHOTOMETRY 

Differential magnitudes for the central star candidate CS1 have been measured using several comparison stars. 
However, the most consistent photometric results and the highest signal to noise ratio was obtained using 
comparison star number 2 for BOS and comparison star number 1 for FTS and LT. Tables 1.1-1.5 list the final 
differential magnitudes for star CS1. 

2 RADIAL VELOCITIES 

Radial velocities (RVs) of the central star candidate CS1 have been measured as relative RVs using cross¬ 
correlation technique and the first spectrum as a template spectrum in each of the two datasets, INT and 
AAT. Relative RVs were then scaled into absolute values by measuring the wavelength of Mgn 4481 A line in 
every spectrum in both datasets. Tables 2.1 and 2.2 list the absolute heliocentric RVs from the INT and AAT 
datasets, respectively. 

3 SYNTHETIC SPECTRA 

Normalized spectra of the central star candidate CS1 have been visually compared to synthetic spectra with 
dispersion 1 A/pixel, computed by Munari et al. (2005) with the Kurutz code. Comparison between the 
measured spectrum and synthetic spectra of various rotational velocities is shown in Fig. 3.1, comparison to 
synthetic spectra with various surface gravities in Fig. 3.2 and various effective temperatures in Figs. 3.3 and 
3.4. Table 3.1 demonstrates the influence of the effective temperature on the ratio between Hei (4471 A) and 
Mg II (4481 A) lines. The effective temperature was determined to vary within a 2000 K range primarily due 
to the measured variable ratio between Hei and Mgn lines (see Table 3.2). The metallicity was determined to 
be larger than the largest available value of metallicity in the library of synthetic spectra, [M/H] > 0.5. The 
comparison between measured and synthetic spectra with the final selection of atmospheric parameters in the 
entire wavelength range is shown in Fig. 3.5. The atmospheric parameters, used for the computation of the 
best-fit spectra, are listed in Table 3.3 and correspond to B9V-B7V (Lang, 1992). 

4 STANDARD SPECTRA 

The measured spectra of the central star candidate CS1 have been compared to a number of standard stellar 
spectra of various spectral types, provided by Bagnulo et al. (2003). The comparison was used to verify the 
~B8V spectral type of CS1, determined by comparisons with synthetic spectra. Fig. 4.1 shows the comparison 
between the measured INT spectrum of the central star candidate CS1 and the selected sample of standard 
spectra of spectral types A7V, B8V and BOV. 
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Table 1.1: Log of differential magnitudes for the central star candidate CS1 using the BOS, 
Liter B and comparison star number 2. 


JD - 2450000 

A B 

mag(2) - nrag(CSl) 

JD - 2450000 

A B 

mag(2) - nrag(CSl) 

5739.816 

0.8217i0.0066 

5751.812 

0.9764i0.0063 

5743.837 

0.8940±0.0058 

5761.734 

1.2197i0.0065 

5743.902 

0.9102i0.0061 

5762.762 

1.1494i0.0061 

5743.938 

0.9067i0.0064 

5764.733 

1.1661i0.0058 

5748.897 

0.8611i0.0064 

5764.741 

1.1679i0.0058 

5748.904 

0.8462i0.0066 

5765.764 

1.2898i0.0057 

5749.768 

0.9490i0.0065 

5768.712 

1.4281i0.0062 

5749.842 

0.9660i0.0062 

5769.716 

1.3524i0.0058 

5749.856 

0.9776i0.0061 

5770.721 

1.4605i0.0064 

5751.798 

5751.805 

0.9816i0.0069 

0.9802i0.0064 

5771.743 

1.4890i0.0060 


Table 1.2: Log of differential magnitudes for CS1 using the FTS, Liter V and comparison 
star number 1. 


JD - 2450000 

A V 

mag(l) - mag(CSl) 

JD - 2450000 

A V 

mag(l) - mag(CSl) 

5838.890 

-2.2670i0.0095 

5997.249 

-1.9800i0.0085 

5842.947 

-2.0290i0.0095 

5998.292 

-2.1220i0.0095 

5843.896 

5844.897 

-1.9250i0.0250 

-1.8890i0.0085 

6010.305 

-2.5360i0.0114 


Table 1.3: Log of differential magnitudes for CS1 using the FTS, Liter B and comparison 
star number 1. 


JD - 2450000 

A B 

mag(l) - mag(CSl) 

JD - 2450000 

A B 

mag(l) - mag(CSl) 

5838.889 

-2.3440i0.0155 

5997.248 

-2.0730i0.0168 

5842.946 

-2.0760i0.0168 

5998.292 

-2.1870i0.0228 

5843.896 

5844.897 

-1.9720i0.0323 

-1.9510i0.0139 

6010.304 

-2.5940i0.0226 


Table 1.4: Log of differential magnitudes for CS1 using the LT, Liter V and comparison star 
number 1. 


JD - 2450000 

A V 

mag(l) - nrag(CSl) 

JD - 2450000 

A V 

mag(l) - mag(CSl) 

6018.721 

-2.5310i0.0116 

6092.672 

-2.5465i0.0113 

6046.594 

-2.0740i0.0095 

6093.683 

-2.5467i0.0118 

6049.576 

-1.8556i0.0080 

6095.463 

-2.5685i0.0151 

6052.568 

-2.1912i0.0109 

6095.562 

-2.5351i0.0113 

6056.665 

-2.3214i0.0103 

6103.648 

-2.3478i0.0115 

6056.668 

-2.3117i0.0102 

6106.541 

-2.4902i0.0144 

6060.608 

-2.1531i0.0268 

6108.567 

-2.5055i0.0138 

6074.515 

-2.3415i0.0118 

6109.574 

-2.4398i0.0121 

6076.508 

-2.4097i0.0113 

6110.522 

-2.2982i0.0111 

6076.528 

-2.3913i0.0111 

6110.525 

-2.3294i0.0120 

6079.493 

-2.4715i0.0127 

6114.524 

-2.0552i0.0093 

6082.700 

-2.5325i0.0111 

6116.584 

-1.9536i0.0091 

6083.600 

-2.5409i0.0113 

6117.558 

-1.9007i0.0092 

6084.640 

-2.5048i0.0111 

6121.565 

-1.8546i0.0091 

6087.645 

-2.5043i0.0111 

6124.565 

-1.9821i0.0117 

6090.694 

-2.5145i0.0113 

6128.497 

-1.7330i0.0100 

6090.697 

-2.5012i0.0110 

6128.500 

-1.7287i0.0093 

6091.668 

-2.5048i0.0105 

6129.575 

-1.8577i0.0127 
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Table 1.5: Log of differential magnitudes for CS1 using the LT, filter U and comparison star 
number 1. 


JD - 2450000 

A U 

mag(l) - mag(CSl) 

JD - 2450000 

AU 

mag(l) - mag(CSl) 

6018.722 

-2.9742i0.0277 

6091.668 

-2.8498i0.0256 

6046.595 

-2.3111A0.0267 

6092.673 

-2.9139i0.0250 

6049.576 

-1.9648i0.0186 

6093.683 

-2.9632i0.0293 

6052.569 

-2.2126A0.0228 

6095.562 

-3.0038i0.0279 

6056.666 

-2.5361A0.0212 

6103.649 

-2.7321i0.0279 

6056.668 

-2.4795i0.0203 

6106.541 

-2.7709i0.0301 

6060.612 

-2.4066i0.0443 

6108.568 

-2.8010i0.0278 

6074.516 

-2.6633i0.0266 

6109.574 

-2.6252i0.0246 

6076.509 

-2.7534i0.0274 

6110.523 

-2.4069i0.0235 

6076.529 

-2.7148i0.0252 

6110.525 

-2.4842i0.0229 

6079.494 

-2.8251i0.0320 

6114.525 

-2.2627i0.0212 

6082.700 

-2.8596i0.0273 

6116.585 

-2.1451i0.0184 

6083.600 

-2.9087i0.0258 

6117.559 

-2.0160i0.0175 

6084.640 

-2.8747i0.0253 

6121.566 

-2.0329i0.0174 

6087.645 

-2.8669A0.0239 

6128.498 

-2.0713i0.0190 

6090.695 

-2.9846i0.0264 

6128.500 

-2.0572i0.0192 

6090.697 

-2.9328A0.0257 

6129.575 

-2.1450i0.0258 


Table 2.1: Log of heliocentric RVs for the central star candidate CS1 using the INT. The 
Erst error bar corresponds to the error of absolute RV scaling by Mg II line, the second error 
bar is the cross-correlation error. 


JD - 2450000 

RV (km s' 1 ) 

JD - 2450000 

RV (km s' 1 ) 

275.4937 

1.48il.31 

676.5633 

29.78ii.31i3.16 

308.4805 

44.58il.31i3.18 

683.5339 

26.17il.31il.65 

320.5166 

-l.03il.31il.90 

704.3787 

30.56il.31i2.39 

354.4195 

38.79il.31i2.25 

731.3648 

17.76il.31i2.66 

384.3285 

19.58il.31il.66 

734.3913 

26.49il.31i2.34 

561.6797 

22.59il.31i5.20 

889.6922 

25.11il.31i2.62 

562.6892 

22.15il.31i5.88 

902.7299 

7.80il.31i4.82 

578.6710 

62.78il.31i2.66 

952.6966 

26.23il.31i3.23 

616.6297 

45.33il.31i2.52 

967.6333 

-5.43il.31i2.59 

631.6013 

9.71il.31i2.43 

972.5833 

-3.34il.31i2.51 

633.5566 

37.80il.31i3.16 

1066.4821 

45.52il.31i2.14 

653.5404 

32.79il.31i2.39 

1303.7116 

55.06il.31i3.81 

675.5400 

39.72il.31i2.40 




Table 2.2: Log of heliocentric RVs for CS1 using the AAT. The first error bar corresponds 
to the error of absolute RV scaling by Mgn line, the second error bar is the cross-correlation 
error. 


.JD - 2450000 

RV (km s" 1 ) 

JD - 2450000 

RV (km s" 1 ) 

1686.14051 

1686.16222 

1686.17858 

1686.21750 

1686.23923 

12.09il.01 

11.46il.01i0.19 

12.13il.01i0.20 

11.32il.01i0.20 

11.14il.01i0.21 

1686.26094 

1686.28265 

1686.30438 

1686.32536 

11.76il.01i0.22 

11.55il.01i0.21 

11.28il.01i0.21 

12.09il.01i0.30 
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Wavelength (Angstrom) 


Figure 3.1: Determination of the rotational velocity. 
The first INT spectrum, obtained on 1996 July 10 
(black), is compared to synthetic spectra with stellar 
rotational velocities ranging from 100 km s" 1 (green) 
to 300 km s" 1 (blue). The basal widths of the ab¬ 
sorption lines were best reproduced by a rotational 
velocity of 200 km s _1 (red). Other atmospheric pa¬ 
rameters are the same for all synthetic spectra: effec¬ 
tive temperature 11000 K, log(g) 4.0 cm s" 2 and solar 
metallicity. The two plotted spectral features are He I 
(4471 A) and Mgn (4481 A). 



Figure 3.2: Determination of the surface gravity. The 
Hd profile is fitted by synthetic spectra with 3.5 
(green), 4.0 (red) and 4.5 cm s" 2 (blue) logarithmic 
stellar surface gravitational acceleration. The best 
fit of the Balmer wings was achieved with log (g) = 
4.0 cm s' 2 (red). The central parts of the Balmer line 
profiles were affected by the emission components and 
were excluded from the fitting process. 



Wavelength (Angstrom) 


Figure 3.3: The effect of the effective temperature on 
the Hei/Mgn ratio for the temperatures given in Ta¬ 
ble 3.1, from 10000 K (green) to 15000 K (blue). A 
higher effective temperature increases the strength of 
the He I line (4471 A), whereas the strength of the 
Mg II line (4481 A) is slightly reduced. The mea¬ 
sured He i/Mg II ratio from the first INT spectrum ob¬ 
tained on 1996 July 10 is best reproduced by the syn¬ 
thetic spectrum with effective temperature 11000 K 
(red line). Note how the Hei line becomes affected 
by the Tin line (4468.5 A) at lower effective temper¬ 
atures. 



Wavelength (Angstrom) 


Figure 3.4: The impact of the effective temperature 
on the Balmer lines, from 10000 K (green) to 15000 K 
(blue). The absorption wings of the Hd profile are well 
fitted by the 11000 K synthetic spectrum (red). This 
confirms the choice of 11000 K, which was made based 
on the Hei/Mgn ratio (see Table 3.2 and Fig. 3.3). 
As for the gravity htting process, the central region 
of the Balmer prohles have not been fitted, because 
of the presence of the emission and self-absorption 
components. 
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Wavelength (Angstrom) 


Figure 3.5: The best fit synthetic spectrum (red) is shown against the measured INT spectrum from 1996 
July 10 (black). The atmospheric parameters of the synthetic spectrum are T e ff=11000 K, log(g)=4.0 cm s" 2 , 
u ro t=200 km s" 1 and [M/H]=0. The discrepancies between the measured and synthetic metal line strengths are 
a result of the high metallicity of the observed star. 


Table 3.1: Hei (4471 A) to Mgn (4481 A) ratio de¬ 
pendence on the effective temperature. The ratios 
are calculated from synthetic spectra with parameters 
log(g) = 4-0 cm s' 2 , [M/H]=0 and u rot = 200 km s' 1 . 


Table 3.2: Estimated effective temperature. The mea¬ 
sured ratio between Hei (4471 A) and Mgn (4481 A) 
and the corresponding effective temperature (see Ta¬ 
ble 3.1) is provided for various INT spectra. The INT 
spectrum, obtained on 1996 October 27, exhibits the 
smallest line ratio, while the largest ratio can be seen 
in the spectrum, taken on 1998 March 16. 


T eS [K] 

He i/Mg ii 

Teff [K] 

Hei/Mg ii 

observation date 

He i/Mg ii 

T eS [K] 

10000 

0.43 

12000 

0.82 

1996 Jul 10 

0.55 

11000 

10500 

0.45 

13000 

1.07 

1996 Oct 27 (min) 

0.51 

11000 

11000 

0.53 

14000 

1.66 

1998 Mar 16 (max) 

1.12 

13000 

11500 

0.65 

15000 

2.70 

average 

0.80 

12000 


Table 3.3: The best ht atmospheric parameters of CS1. The effective temperature, stellar 
surface gravity and projected rotational velocity have been fitted with the synthetic spectra. 
Metallicity could not be fitted, because it exceeded the values, provided in the library of 
synthetic spectra. 


T eS [K] 

11000-13000 

log(g) [cm s' 2 ] 

4.0 ±0.25 

u rot • sin( j) [km s' 1 ] 

200^25 

[M/H] 

> 0 
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Wavelength (Angstrom) 


Figure 4.1: Comparison between our measured INT spectrum of CS1 (black) and standard spectra A7V (top), 
B8V (middle) and B6V (bottom) shown in red. A7V: Note the absence of Sin lines (3856 and 3863 A), Hei 
(4026 and 4471 A) and presence of Mm (4033 A), Cal (4227 A) and Tin (4468 A). The Balmer prohles exhibit 
much stronger gravitational broadening than in our measured spectra. B8V: Sin and Hei lines are present, 
whereas Ti II (4468 A) is not because of the higher temperatures of B8 stars. The measured depth ratio between 
Hei (4471 A) and Mgn line (4481 A) is perfectly matched. Also note the comparable amount of gravitational 
broadening in the Balmer absorption lines. The narrowness of the other lines may be disregarded because of the 
smaller projected rotational velocity of the B8V standard star. B6V: The 2000 K higher effective temperature, 
compared to B8, results in highly pronounced Hei lines (4026 and 4471 A), whereas the lines of the heavier 
elements are no longer present. The larger radius and the corresponding smaller gravitational broadening of a 
B6V star mean that the absorption wings of the Balmer prohles are not as well matched as for the B8V star. 
The most noticeable discrepancy, the large ratio between the He I and Mg II line, suggests that the spectral type 
of the observed star has to be later than B6. 
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